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Abstract 

Earlier we introduced (M.I. Trofimov, E.A.Smolenskii, Application of 
the Electronegativity Indices of Organic Molecules to Tasks of Chemical 
Informatics, Russian Chemical Bulletin 54(2005), 2235-2246. 

|http : //dx . doi . org/10 . 1007/slll72- 006-0105-6) 

effective recursive algorithm for graph isomorphism testing. In this 
paper we describe used approach and iterative modification of this algo- 
rithm, which modification has polynomial time complexity for all graphs. 
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1 Introduction. 

As 30 years ago [4j [9] , in spite of efforts of many researchers, whether the 
problem of testing of graph isomorphism is NP-complete is still open. There 
is an interesting explanation of this fact [3j. The authors noted that proofs 
of NP-completeness seem to require a certain amount of redundancy, which 
redundancy graph isomorphism problem lacks. For example, in the case of 
subgraph isomorphism search the same result may be observed even if some 
edges of given graph will be deleted or some new edges will be added. In 
contrast, if any edge will be added to one of two isomorphic graphs (or if any 
edge will be deleted), then the graphs will no longer be isomorphic. So the graph 
isomorphism problem is not typical NP-complete problem. At present many 
effective algorithms for graph isomorphism testing were proposed and there are 
number of investigations with attempts to prove polynomial complexity of some 
of these algorithms. For example, we note preprint [2]. 

This article deals with an analysis of complexity of iterative modification of 
graph isomorphism test algorithm, which had been described earlier The 
initial recursive algorithm was developed for effective solution of number of hard 
tasks of computer chemistry, however, here we will not consider the chemical 
applications of this approach, we only note that in contrast with many other 
formal approaches, which are used for graph theory applications in chemistry 
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today, this approach is based on fundamental physico-chemical notions about 
structure and properties of a substance. 

To describe the algorithm we will use Pascal-like pseudo-code, for example, 
similar pseudo-code had been used in well-known classical books [TJ [T3] . 



2 The general principles of the approach. 

Further, without loss of generality of the task, we will consider undirected con- 
nected graphs without loops [SJ [5]. For the algorithm we will use weights of 
vertices, which weights may be computed from the following equation |llj : 



±-^ Xj + bi, (1) 



where di is degree of vertex i; 

j is vertex number for all i-adjacent vertices (i.e. vertex j is i-adjacent if an 
edge exists in the graph); 

bi is absolute term, positive number (the cases, where 6^=0 will be noted 
specially) . 

I.e. the weight of every vertex i is arithmetic mean of i-adjacent vertices 
weights and initial weight of ith vertex Bi = (di + l)bi. Having written similar 
equations (UJ for all n vertices of the graph, we get the system of n linear equa- 
tions, with solutions xi,x%, ...,x n . Every system ([1]) always has single solution. 

In matrix form system ([TJ) can be written as 

A-X = b, 

where A = E — L; 

E is the unity matrix of corresponding size; 
M = 1 1 my 1 1 is adjacency matrix of the graph; 
L = ||* y || = JV-M; 



ifmij = 1, 
otherwise; 

I otherwise. 

Consider the matrix A~~ x = Since lim L h = 0, we obtain Neumann 

h— >oo 



series: 



A' 1 = (E- L)- 1 = J2 Lk = E + J2 Lk - 
From that follows: 



h=0 h=l 
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l<kij < 2 if ' i = j, 
< kij < 1 otherwise. 



(2) 



Two graphs G and G with the adjacency matrices M and M', respectively, 
are isomorphic if and only if there exists the permutation matrix 



P = \\Pi 



where = 



1 if * J, 
otherwise; 



such that ([5], p. 151): 

M' = P- 1 - M-P (3) 

Note following properties of P: 

1) P is an orthogonal matrix, i.e. P _1 = P T ; 

2) if G = G" and i — > j, where i is vertex of G, j is vertex of G', then 
~e~i = P • ej and En = P ■ Ejj ■ P -1 , where all coordinates of the vector 
are equal to except for the ith coordinate equal to 1 and the element on the 
crossing of the ith row and ith column of the matrix En is equal to 1, while all 
other elements of En are equal to zero. 

Let's define di,d[ vectors with the coordinates di : d' i: respectively. In the 
matrix form one can write d = M ■ e and d' = M' ■ e, where e is unity vector. 
From that follows d' = P^ 1 ■ d. Let's define D, D' as the diagonal matrices with 
the diagonal elements equal to di : d' i: respectively. Then iV = (E + D)^ 1 and 
N' = (E + D') -1 , where E is the unity matrix of corresponding size. Also D 
and D' matrices can be expressed through the di and d. t vectors: 

D = Y / Eu-d-ef,D' = Y / E n -d'-ej. 

i j 

Consider two systems of linear equations ((1) for G and G', respectively: A-X = b 
and A' ■ X~' = t . 

Assertion 1. Let G = G' and let P be the permutation matrix. Then if 
t = P- 1 ■ b, then X' = P- 1 ■ X. 

Proof. Let's find a relation between the matrices A and A': 

P- 1 ■ D ■ P = P" 1 • (£ En ■d-ef)-P = J2P- 1 -E u -d-ef • P = E P" 1 • 

E u ■ P ■ P- 1 ■ d ■ (P- 1 ■ e 4 ) T = £ E n ■ d> -ej=D>. 

j 

From that follows: 

N' = (E + D)- 1 = (P- 1 ■ E ■ P + P- 1 ■ D ■ P)- 1 = (P- 1 ■ {E + D) -P)- 1 = 
P- 1 ■ (E + D)- 1 ■ P = P- 1 ■ N ■ P 
and also: 

A' = E-N'-M' = P- l -E-P-P- l -N-P-P~ l -M-P = P- 1 -{E-N-M)-P = 
P- 1 ■ A- P, A'- 1 = P- 1 ■ A- 1 ■ P. 
Finally: 
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X' = A'- 1 ■ b' = P- 1 ■ A- 1 ■ P- P- 1 -b = P- 1 ■A- 1 -b = P- 1 ■ X, 
which was to be proved. ■ 

Assertion 2. For any two graphs G and G', which graphs have n vertices, 
and for any permutation P £ S n , if for all b, and b = P^ 1 ■ b , and for 

solutions X and X of system (1) we have X = P^ 1 ■ X, then P is the graphs 
isomorphism. 

Proof. Let M be the adjacency matrix of the graph G , let M' be the 
adjacency matrix of the graph G" , and M" be the adjacency matrix of the 
graph G" . Let M" = P^ 1 ■ M ■ P and b = b . Then according to definition 
of graph isomorphism, from M" = P _1 ■ M ■ P we have G = G" . So, according 
to Assertion 1 from b = P^ 1 ■ b we have X = P^ 1 ■ X. Taking into account 
X = P~ x ■ X, we get X — X . Also A ■ X = b, so from b = b we have 
A" ■ X" = A' ■ X'. Taking into account it = ~x" , we have A' = A". This 
means that M' = M". Taking into account M" = P^ 1 ■ M ■ P, finally we have 
M' = P- 1 - M-P and G = G' . ■ 

It is self-evident that n linear independent vectors b are sufficient to corre- 
spond to Assertion 2. From that we have following assertion: 

Assertion 3. For any two graphs G and G' , which graphs have n vertices, 
and for any permutation P € S n , if for n linear independent vectors b and for 
n vectors b = P^ 1 ■ b , and for solutions X and X of system (1) we have 
X = P^ 1 • X , then P is the graphs isomorphism. 

According to Harary's definition: "Two points [vertices - MT] u and v of the 
graph G are similar if for some automorphism a of G, a(u) = v." [S], p. 171. 
We expand this definition for vertex u of the graph G and vertex v of the graph 
G'\ 

Definition 1. If G = G' and v ■<-> u is possible mapping, then u and v are 
similar M 

3 The algorithm. 

Let graphs G, G' are isomorphic. For selected pair of vertices (i, j),i G V, j G V 
with the same degree we build bigraph Hij — (V, V, C/y), where V is set of 
vertices of graph G; V is set of vertices of graph G'; C/y is set of edges of 
bigraph Hij, that vertex p £ V connected with vertex q 6 V, if deg(p) = deg(g) 
and kip — k'j q . The number of edges in is not more than n 2 . For vectors 

Sj, e.j solutions of systems (p} : Xi\p] = XAq], if there is edge (p, q) £ J7y; and 
Xi[p] ^ if there is not edge (p,q) £ JJy. Due to @, every vertex from 

pair in bigraph ify is incident to only one edge (the edge is edge 
From ([2]) we also can see that all fcy ^ 0. 

Let function T(U) returns a transversal (subset of edges of bigraph H = 
(V, V', U)), if such transversal exists (if there are a few transversals, then the 
function returns any transversal), and empty set otherwise. 

Now we can consider the following procedure: 
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procedure Pl(i,j, R); 

1) For input pair of vertices i £ V, j £ V we build bigraph H^, as 
described above. 

2) For every edge (p, q) from Uij (except edge we build bigraph H pq , 
as described above. If T(Uij f] U pq ) ^ 0, then replace Uij :— Uj f] U pq , else 
remove edge (p, g) from C/y . 

3) If T(Uij) 7^ 0, then found transversal gives isomorphic map for given 
graphs (the result map will be returned via R)M 

Assertion 4. If input vertices are similar (G = G' ) then procedure PI finds 
isomorphism. 

Proof. 1) Procedure PI tries to find isomorphism via one-to-one corre- 
spondence with vertex weights of given graphs. These weights are defined by 
functions of initial weights of vertices and so the functions have to be bijective 
(to exclude degenerate cases), i.e. the necessary condition is one-to-one corre- 
spondence with weight of every vertex p and initial weight of every vertex q. 
Let's prove that used weights correspond to this condition. The weight of vertex 
p is 

n 

t=i 

Suppose that 

C q = kpth = const, 

tjiq 

then 

Xp = f{bq^j = kpqbq + Cq. 

This function is bijective when k pq ^ 0. As was noted above, from ([2]) we 
get that all k pq ^ 0. Hence the result. 

2) Let's prove that if in the result of PI T(Uij) ^ 0, then this transversal 
is single. As was noted above, every vertex from pair (p, q) in bigraph H pq is 
incident to only one edge and this edge is edge (p, q). So, for vertices p, q there 
is only one edge (p, q) in intersection Uij f] U pq , and so after edges replacement 
Uij := Uj P| Upq in bigraph Hij: vertex p is incident to only one edge (p, q) and 
vertex q is incident to only one edge (p,q). After the loop for all edges (p, q) 
from Uj (step 2 of PI) we get bigraph Hij, where every vertex is incident to 
only one edge, so transversal T(Uj) is single. 

3) Let's prove that if in the result of PI T(Uj) ^ 0, then this transversal 
maps isomorphism of given graphs. As was noted above: Xi[p] — Xj[q], if 
there is edge (p, q) in bigraph Hij. In this case T(Uj) is common transversal 
for Uj and for U pqi i.e. for every edge [u, v) £ Uij we can find edge (u, v) £ 
Upq. Consider permutation P, where p — > q. if edge (p,q) exists. For n linear 
independent vectors e%, ...,e„: let b p — e p , p = 1, ...,n; b = P^ 1 ■ b. Then we 
get X = P^ 1 ■ X. According to Assertion 3: P is isomorphism. 
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4) Let in loop 2 of PI the first n edges of bigraph Hij give isomorphism 
(pi : i <H- j, r\ <->■ r[, ...,r n _i f-> r^^. In this case, as was proved (section 2 of 
this proof): every vertex is incident to only one edge. 

Let next step of the loop finds edge (p,q), where vertices p, q are not sim- 
ilar, and k lp = k' jq . Then T(Uijf)U pq ) = 0. Otherwise T(Uijf]U pq ) would 
be common transversal for U rir > , U rn _ irn _ 1} and any isomorphism p O q 
would be possible, and so vertices p, q would be similar, which contradicts our 
assumption. 

Let next step of the loop finds edge (p,q), where vertices p, q are similar. 
Then if T(Uij |") U pq ) ^ 0, then another isomorphism <p2 will be found, otherwise 
the edge (p, q) will be removed and initial isomorphism (f>\ will be saved in result 
of procedure PI. 

Taking into account commutative and associative properties of intersection 
we have in conclusion that any isomorphism will be found for every search order 
of loop 2.1 

Now we can consider the following algorithm, where function Verify tests 
found mapping by procedure PI. The mapping may be tested via equation ([3]): 
if the mapping is isomorphism, then the function returns true, otherwise the 
function returns false. 



Algorithm 1. 



Input: graphs G, G' 



if (numbers of vertices of G and G'are different) 

or (numbers of edges of G and G'are different) then 
goto 1; 

Calculate A' 1 , A'^ 1 ; 

for j:=l to n do 
begin 

Pi(i,i,#); 

if Verify (G,G',R) then 
begin 

write('G = G' '); 
goto 2; 

end 

end; 

1: write ('G £ G"); 
2: ;■ 



Assertion 5. Algorithm 1 finds isomorphism for given graphs. 

Proof. If graphs G, G' are isomorphic, then there is vertex j (j = 1, n) 
in graph G', which vertex j is similar with vertex 1 of graph G . In this case 
according to Assertion 4 procedure PI finds isomorphism for input pair (l,j) 
and function Verify returns true, the algorithm prints a message "G = G'" and 
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stops. 

If graphs G, G' are not isomorphic, then if numbers of vertices of G and G'are 
different or if numbers of edges of G and G'are different, then the algorithm 
prints a message "G ^ G'" and stops. Otherwise procedure PI may find any 
map R, but function Verify returns false and after n steps of the loop "for j:=l 
to n do" the algorithm prints a message "G p= G'" and stops. ■ 

Proved Assertion 5 suggests that graph isomorphism problem may be re- 
duced to particular case of common representatives system finding. 

4 Complexity of the algorithm. 

The complexity of calculation of matrix A" 1 may be estimated [7] as 0(n 3 ). The 
number of edges in Hij for procedure PI is not more than n 2 (see above), hence 
the complexity of bigraph building may be estimated as 0(n 2 ) and PI uses not 
more n 2 +l bigraphs. So the complexity of PI may be estimated as 0(n 4 ) and 
the main effort of Algorithm 1 is n calls of PI in the worst case (function Verify 
may be implemented much more effective). Hence total complexity is 0(n 5 ). 

5 Conclusion. 

Side by side with general theoretical result of this work we would like to note its 
applied result. Sorted solutions of system (1) are graph invariants (topological 
index). As we can see, the index has high discriminating capability, so it may 
be useful for many practical tasks outside the graph isomorphism problem [TT| , 
sometime with combination with other indices |12| , for example, with modifica- 
tions of Hosoya's index |10j . The properties of this index discussed here may be 
useful for these tasks as well. The researches of system (1) had been started in 

1987 [g. 

Possible corrections and supplemental files will be available via 
http : / /mt2 . comtv . ru/| 
and via 

|http : / / webcenter . ru/~mtrof imov/ 

The password to unzip exe-files is 

hH758-kT402-N3D8a-961fQ-WJL24 
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